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Inelastic cotunneling current and shot noise of an interacting quantum dot with 

ferromagnetic correlations 
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We explore inelastic cotunneling through a strongly Coulomb-blockaded quantum dot attached 
to two ferromagnetic leads in the weak coupling limit using a generic quantum Langevin equation 
approach. We first develop a Bloch-type equation microscopically to describe the cotunneling- 
induced spin relaxation dynamics, and then develop explicit analytical expressions for the local 
magnetization, current, and its fluctuations. On this basis, we predict a novel zero-bias anomaly of 
the differential conductance in the absence of a magnetic field for the anti-parallel configuration, and 
asymmetric peak splitting in a magnetic field. Also, for the same system with large polarization, 
we find a negative zero-frequency differential shot noise in the low positive bias-voltage region. All 
these effects are ascribed to rapid spin-reversal due to underlying spin-flip cotunneling. 
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I. INTRODUCTION 

Recently, there has been extensive investigation of 
spin-dependent tunneling in nanoscale semiconductor de- 
vices coupled to ferromagnetic leadsi. It is well known 
that transport through a quantum dot (QD) includes two 
distinct mechanisms: sequential (first-order) tunneling if 
the QD level is in resonance with the Fermi levels of 
the electrodes, and higher-order tunneling (cotunneling) 
when the QD level is far removed from resonance, with 
sequential tunneling exponentially suppressed. Second- 
order cotunneling in the weak tunneling regime has 
already been experimentally observed at temperatures 
above the Kondo temperature 2 . So far, most studies con- 
cerning spin-polarized tunneling through a QD have been 
focused on the sequential regime^ and on enhanced co- 
tunneling in the strong tunneling regime, i.e. the Kondo- 
type transport^* 6 -. In particular, the theoretical predic- 
tion of Kondo peak-splitting in a quantum dot in the 
presence of parallel (P) spin-polarized leads^ has been 
confirmed in recent experiments 6 . 

In this paper, we systematically analyze cotunneling 
through a strongly Coulomb-blockaded QD attached to 
two ferromagnetic leads in the weak tunneling regime at 
temperatures above the Kondo temperature. We develop 
an explicit expression for QD magnetization, which is be- 
lieved to be of essential importance in understanding the 
nonequilibrium Kondo effect, albeit that the formula is 
derived in second-order perturbation theory 7 ^. Our cal- 
culation of cotunneling current predicts a novel zero-bias 
peak (ZBP) of differential conductance in the absence of 
a magnetic field, and an asymmetric peak splitting in 
a magnetic field for the anti-parallel (AP) configuration. 
Our discussion shows that these interesting features arise 
from spin-flip processes which induce rapid spin-reversal. 
An analogous ZBP in cotunneling was found in an ear- 
lier work by Weymann^. Here our explicit analytic ex- 



pressions show similar results of the peak height and its 
temperature dependence as in Ref . ©. 

Moreover, we present an analytical investigation of 
zero-frequency shot noise in the cotunneling regime, 
which has been previously studied only in the sequen- 
tial regime^. We predict that the differential shot noise 
at low bias-voltage region is heavily dependent on the 
polarization of electrodes in the AP configuration. 

The remaining parts of the paper are arranged as fol- 
lows. In Sec. II, we present the physical model and 
theoretical formulation used in this paper. We give the 
Bloch-type dynamical equations derived from a generic 
quantum Langevin equation approac h 10 ! 11 ' 12 to describe 
the time evolution of the QD spin and all relevant cor- 
relation functions. In Sec. Ill, we discuss the nonequi- 
librium magnetization of the QD in detail and point out 
that the spin-flip cotunneling will induce a rapid spin ac- 
cumulation in the AP configuration. It is our main result 
for theoretical explanation of the novel behavior of the 
differential conductance and zero-frequency shot noise, 
whose explicit expressions are derived based on the lin- 
ear response theory in Sec. IV. Numerical evaluations 
and elaborate discussions are also provided in this sec- 
tion. Finally, a brief conclusions are given in Sec. V. 



II. THEORETICAL MODEL AND 
FORMULATION 



Transport through an interacting QD with a single 
energy level ed connected to two ferromagnetic leads is 
modeled by a single impurity Anderson Hamiltonian. To 
describe cotunneling through the QD in the strongly 
Coulomb-blockaded regime (the on-site Coulomb inter- 
action U — > oo), in which the negative with <C \ii v \ 
(fi n is the chemical potential in lead n and the symmet- 
rically applied bias voltage V has hl = —[ar = eV/2) is 
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satisfied, so that charge fluctuations are completely sup- 
pressed, we perform a mapping of the Anderson Hamil- 
tonian onto the subspace with one electron in the dot by 
a Schrieffer- Wolff unitary transformation^, i.e. onto the 
effective Kondo Hamiltonia n 12 ' 14 : 



H —H + Hi, 



(i) 
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where c^ k(r (c^ko-) is the creation (annihilation) operator 
for electrons with momentum k, spin-er and energy e^k 
in lead r\ (= L, R), S z ^ are Pauli spin operators of elec- 
trons in the QD, and J n ri' is the Kondo exchange coupling 
constant, -ffdir is the potential scattering term, which is 
decoupled from the electron spin due to the number of 
electrons in the dot level being one. As a result, this term 
has no influence on the dynamical evolution of the elec- 
tron spin and behaves only as a direct bridge to connect 
the left and right leads and consequently to contribute 
a current being independent of the dynamics of the QD. 
For an Anderson model with symmetrical coupling to the 
leads t, we have J vn > = 2J = t 2 / £ rf- 

In this paper, we consider the ferromagnetism of the 
leads by means of a spin-dependent, flat density of states, 
P m - PL] = PR] = (l+p)p ; Pl\ = PR{ = (1 -p)Po 
for the P configuration, and pl] = Pri = (1 + p)po', 
Pli = p R ^ = (1 —p)po for the AP configuration, with the 
degree of spin polarization p (\p\ < 1) for both leads. 

The spin-splitting, A, of the QD involves two contri- 
butions: an ambient magnetic- field £>-induced Zeeman 
term, Az = gpB-B and an effective splitting, A P (ap) = 
D J2 V J vv(Pv\~Pvl) = 2pD(J LL ±J RR )p Q , caused by fer- 
romagnetic correlations of the two leads (D is the band 
width of the leads). Based on this Hamiltonian, we find 
that (1) the differential conductance, dl/dV of the co- 
tunneling current exhibits double peaks with width 2A p 
as a function of voltage, V, for P alignment of the lead 
magnetizations even without an ambient magnetic field; 
(2) the peak-splitting in the dl/dV vs. V curve with- 
out a magnetic field is also predicted for the AP con- 
figuration only with asymmetric couplings Jll ^ Jrr, 
but it has a largely reduced width, 2 Aap , in comparison 
with the P case; (3) in both cases the splitting can be 
removed by properly tuning the direction and strength 
of the ambient magnetic field, to induce vanishing of the 
total spin-splitting A. These results provide a good qual- 
itative explanation of the experimental findings in Ref . 0, 
and we can further clarify this issue in that the small but 
nonzero splitting found in the asymmetric AP case also 
stems from the lead-magnetization-induced spin-splitting 
in the QD. 



We employ a generic Langevin equation approach to 
analyze this system in the weak cotunneling limit 1 - i 11 ' 12 . 
In our derivation, operators of the QD spin and the reser- 
voirs are first expressed formally by integration of their 
Heisenberg equations of motion (EOM) , exactly to all or- 
ders of J m ' . Next, under the assumption that the time 
scale of decay processes is much slower than that of free 
evolution, we replace the time-dependent operators in- 
volved in the integrals of these EOM's approximately in 
terms of their free evolutions. Thirdly, these EOM's are 
expanded in powers of J vv > up to second order. To this 
end, a Bloch-type dynamical equation is established to 
describe the time evolution of the QD spin variable as, 



S' = -4C P(A P)(A) S z + 2i? P(AP) (A), 



in which 
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Rp(u) = - {g LL + g RR + 2g LR ) (1 - p )u, 
for the P configuration, and 
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+Kg LR (l+p 2 )uj + 2irg LR pV, 
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,(3a) 
(3b) 

;9lrT 

,(4a) 
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for the AP configuration, with g vv > = J^Pq and ip(x) = 
x coth(ai/2). Once again, we observe that the direct tun- 
neling term has no contribution to dissipation of the 
QD spin variables. Throughout, we will use units with 
H = ks = e = 1. 

Based on the dynamical equation, Eq. @, we iden- 
tify the cotunneling-induced spin relaxation rate as 

Pf API 

\jT x = 4Cp(ap)(A), which stems completely from 
spin-flip events involving both single-barrier processes 
[Fig. 3(c)-(h) in Ref. [l2| and double-barrier processes 
[i.e. electron-transferring (ET) processes] [Fig. 3(m)-(p) 
in Ref. [l2j]. For the case of fully spin-polarized leads 
{p = ±1) in the P alignment, cotunneling is unable to flip 
the electron spin in the QD, thus implying that there is 
no spin relaxation, i.e. 1/T\ P = [C(R)p = at p = ±1 
from Eq. I]3a[)]. For this case, we also have S z = 0, mean- 
ing that the QD spin orientation remains unchanged from 
its initial state. In contrast to this, in the case of fully AP 
spin-polarized leads, spin relaxation rate is nonzero and 
it arises solely from ET events [Fig. 1(c) and (d) below]. 



III. NONEQUILIBRIUM MAGNETIZATION 

The nonequilibrium local magnetization, Sp^pj , of the 
QD in P(AP) configuration is readily obtained using the 
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steady-state solution of Eq. ([2]) as 



S£°°(A,V) 



Si? (A, V) 



J?p(A) = 

2C P (A) . . /A 

(9ll + 9rr) ply 
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The P magnetization formula, Eq. (O, is identical with 
previous theoretical results for normal leads^&i^, and it 
is worth noting that the polarization of the reservoirs is 
fully reflected in the effective spin-splitting, Ap. If an 
ambient magnetic field is applied to induce vanishing of 
the spin-splitting, A = 0, we have S§,°°(0, V) = 0. 

However, the local magnetization in the AP configura- 
tion is significantly different. Our calculated results are 
shown as functions of bias- voltage in Fig. 1. The rele- 
vant parameters in our calculations are: qll = gRR = 
9lr — 9 — 0.25 x 10~ 2 ; D = 10 (corresponding to Kondo 
temperature T K = 4.5 x 10~ 5 D and A P = Q.2pD), and 
the temperature is T/D = 0.002. In the absence of an 
external magnetic field, S%p is a two-valued function, 
±p/{l + p 2 ) if \V\/T > 1, and zero net spin of the QD 
occurs only at V — when p ^ [Fig. 1(a)]. Nonzero lo- 
cal spin is generated by the application of a bias- voltage, 
indicating spin accumulation in the QD under nonequi- 
librium conditions, and its orientation is rapidly flipped 
when the direction of the bias-voltage is reversed. This 
interesting behavior can be physically explained by spin- 
flip elastic cotunneling processes, as shown in Figs. 1(c) 
and (d). If V > [Fig. 1(c)], the dominant cotunneling 
process is: first a spin-down electron in the QD tunnels 
out to the right lead, then a spin-up majority electron 
successively enters into the QD. This process continually 
pumps spin-down electrons out of the QD, while injecting 
spin-up electrons into the QD, giving rise to an up-spin 
pilcup, until a steady-state is reached. On the contrary, 
for V < 0, the reverse process takes place, causing a 
down-spin accumulation [Fig. 1(d)]. Another interesting 
consequence of this spin-flip cotunneling is the appear- 
ance of a sharp peak at V = in the vs - V curves, 
which leads to a novel ZBP in differential conductance 
for the AP configuration, as shown in Fig. 2(b). Simple 
calculation shows that the peak height of at V = 
is proportional to the polarization p of electrodes, but is 
inversely proportional to the temperature T: 



While for \V\/T oo, we have = 0. These two 

relations are useful in analyzing ZBP height of differential 
conductance. The temperature behavior of is shown 
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FIG. 1: Nonequilibrium QD magnetization, S 12 , vs. bias- 
voltage, V, for the AP configuration at (a) Az = 0, and (b) 



Az^O. Insets: 
cotunneling processes. 



dV 



V curves, (c) and (d) ET spin-flip 



AT' 



(7) 



A finite ambient magnetic field splits the peak in 
into two peaks at V = ±A with unequal heights 
[Fig. 1(b)]. This occurs because spin-flip cotunnel- 
ing becomes an energy-consuming (inelastic) process in 
a nonzero magnetic field, thus requiring a sufficiently 
strong external bias-voltage for activation. According 
to Eq. ((6|), the spin-reversal takes place approximately 
at A/V = —p. Moreover, the application of a positive 
magnetic field makes the cotunneling shown in Fig. 1(c) 
prevail over that in Fig. 1(d). Correspondingly, the peak 
at V — A is suppressed, but the peak at V = — A 
is significantly enhanced if p > 0. The spin reversal 
property discussed above has potential for application in 
quantum computing and quantum information process- 
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ing. Finally, in the absence of a bias voltage, the lo- 
cal magnetization reduces to the equilibrium expression 
(A,0) = ±tanh(A/2T). 



C200 

D P(AP) 



IV. CURRENT AND SHOT NOISE 

The tunneling current operator of the QD is defined 
as the time rate of change of charge density N v = 

Eka c lka c v ka in lead Jr,(t) = N v . From linear- 
response theory we have 



J — oo 



(8) 



where the statistical average (• ■ • }o is performed with re- 
spect to two decoupled subsystems, QD and reservoirs. 

Nonequilibrium quantum shot noise is another active 
subject in mesoscopic physics, because the current cor- 
relation function can provide further information about 
electronic correlations which is not available from con- 
ductance probes alon o 15 > 16 . The noise spectrum is de- 
fined as the Fourier transform of the current-current cor- 
relation function, S nv >(T), and it too can be calculated 
using linear-response theory: 



S, 



</'/' 



dre^h[SJ n (t),SJ n ,(t')} + ) , (9) 



with 5J v {t) = J v {t) - (J v (t)). 

Considering that the interaction Hamiltonian Hi in 
Eq. fl} contains two components, the direct tunneling 
term H^ lr (the direct tunneling channel) and the cotun- 
neling term (the indirect tunneling channel), we can di- 
vide the current and shot noise into three contributions, 



:[(1-P)M^)-(1 + P)M- 



T 



-)],(12a) 



(l-/)2>(^) 

+r[(i - p)V(^) + (l + P )\^±^)) 

+2S5£[(1 - p) 2 (V - A) - (1 +pf{V + A)](12b) 



for the AP configuration. Further elaborate analy- 
ses clarify that the first linearly bias-voltage-dependent 
terms in Eqs. (jlip and (|12p result from spin-conservative 
ET cotunneling processes, whereas spin-flip ET events 
are responsible for the other contributions^ 2 -. It is note- 
worthy that the first terms on the right hand sides of 
Eq. (fTTa|) [ (fT2a|) ] and Eq. (fT2a|) [ (pb) ] obey the nonequi- 
librium fluctuation-dissipation (NFD) relation, while the 
other two terms represent the generalized NFD relation 
due to energy-consumption cotunneling processes^. The 
direct tunneling yields 18 : 



4/1T9LR - ^{l+p 2 )V, 

9lr 
9o 



(13a) 



S F L f(0)/ng LR = 4^(l+p 2 )7VQ, (13b) 
9lr. 1 



for the P configuration, and 



= i—(i-P 2 )v, 

9LR 

9o 



(14a) 



■ V, 



S£ P L d (0)/ng LR = 4-^(1 -p 2 )T<p{-), (14b) 

9LR 1 

for the AP configuration with go — JqPo — 9lrI^- Be- 
cause the direct tunneling is a spin-conservative process, 
it only interfere with the non-spin-flip cotunneling pro- 
cesses. The interference effect thus contributes to the 
current and shot noise only for the P configuration: 



/ = 7°+/ d + / in , (10a) 

s LL (u) = s^h + ^h + s^m, (iob) 

due to the indirect channel, direct channel, and the in- 
terference effect between the two channels, respectively. 
After lengthy but straightforward calculations, we ob- 
tained explicit expressions for steady-state cotunneling 
current and for frequency-independent auto-correlation 
shot noise via the indirect channel: 

4/7rg LR = (l+p 2 )V + 2(l-p 2 )V + 2TS§.°° 
x(l-p 2 )[^(^)-^-t% (11a) 



T 

,2 > 



> L mh9LR = (i+p 2 )M-) + (i-p 2 )t 

+ - 4SjT(l -P 2 )A,(Hb) 

for the P configuration, and 

I%/*9lr = (l- P 2 )V+[(l-p) 2 (V-A) 
+ (1 +p) 2 {V + A)] + 2T5j^? 



9LR 



(15a) 



S F L ?(0)/ng LR = 16 *L pS *°°T<p(-). (15b) 

V 9LR J- 

Setting p = 0, we obtain the cotunneling current formula 
for the normal-lead/QD/normal-lead system [Eq. (36) in 
Ref. [l2|^. and the corresponding frequency-independent 

'zoo] 

p 



shot noise S LL (0) is given by [S zoa = S£°°] 



Sll(P)/*9lr = 2IV(~) - 4S ZOO A 



-T[p 



■ V - A , 
-f—j 



■ V + A 
T 



)].(16) 



di. 



We plot the differential conductance, Gp(ap) = ^ P Jy F) 
(in units of 47r^), and the tunnel magnetoresistance, 
TMR=(G P - G A p)/G A p, vs. bias-voltage in Fig. 2. 
For the P alignment, we find a characteristic jump at 
V = ±A even at zero ambient magnetic field, Az = 0, 
due to the lead polarization-induced spin-splitting, Ap. 
But for spin-polarized leads with larger polarization, 
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FIG. 2: Differential conductances and TMRs as functions o 
bias-voltage under zero (a)-(c) and nonzero (d)-(f) externa 
magnetic fields for various polarizations p. The other param 
eters are the same as in Fig. 1. 



which is exactly the same as the result in Ref. [||. Fur- 
thermore, we examine the temperature dependence of the 
ZBP as shown in Fig. 3(b). It is found that increasing 
temperature will enhance the width of the ZBP, but re- 
duce the peak height a bit, which is qualitatively consis- 
tent with the numerical calculations in Ref. [see their 
Fig. 1(b)]. Interestingly, the novel ZBP in the AP con- 
figuration is robust over a rather wide region of temper- 
ature. It is worth emphasizing that our investigation is 
valid for the cotunneling through a strongly Coulomb- 
blockaded QD (U — > oo) in the weak-tunneling limit, 
thus the contribution of the charge fluctuation effect is 
completely excluded in the present analysis. We believe 
that this is the reason for the difference between our anal- 
vsis and the numerical calculations of Ref. 




o 
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p = 0.8 and p = 1 in Fig. 2(a), the peak-splitting dis 
appears, as can be understood from the fact that th< 
spin- flip process can hardly (p — 0.8) or fully not (p = 1 
occur in these specific cases, thus no additional chan 
nel is opened for ET even when V > A, in compari 
son with smaller polarization p < 0.8^. For the case o 
vanishing spin-splitting A = 0, the linear conductance 
is readily given by 



no 



dip , 

-dv W- 



(%/wg = 4. 



(17 



On the contrary, for the AP configuration, a ZBP 
emerges. Obviously, this behavior stems mathematically 
from the sharp peak in around V = [Fig. 1(a)] 
according to Eq. (|12ap . Consequently, one can conclude 
that the underlying rapid spin-reversal due to spin-flip 
elastic cotunneling is responsible for the novel ZBP from 
a physical point of view. With the help of Eq. |0, we 
obtain the linear conductance, Ga P = d l$? |v_»n. as 



V 



WTp dSlp 



dV 



4(1 -p 2 ). (18) 



We can also calculate the differential conductance Gap 
at the limit of V/T -> oo: 



G% P /ng = 3+p 2 -8pS z AP \ V - 
It is easy to see that G AP > 



4(1 -P 2 ) 
l+p 2 



(19) 



l_r A p 



is always satisfied, 
which reflects the ZBP. Moreover, the relative height x 
of the ZBP is defined as 



/-r0 
Lx A p 



Cr Ap 



G 



(20) 



FIG. 3: Temperature dependences of ^ r (a) and differen- 
tial conductance (b) as functions of bias-voltage in the AP 
configuration with zero magnetic field and p — 0.5. 



In combination with the peak-splitting in the P con- 
figuration, the ZBP in the AP configuration will cause a 
deep dip in TMR at V — 0, as shown in Fig. 2(c). Analo- 
gously, an ambient magnetic field leads to peak-splitting 
since inelastic spin-flip cotunneling requires sufficient en- 
ergy for activation. Contrary to the P alignment case, the 
two peaks in have unequal heights due in part to 



the asymmetry in with nonzero Az [Fig. 1(b)]. Cor- 
respondingly, the magnetic field changes the behavior of 
TMR [Fig. 2(f)]. 

Richly detailed structures are observed in the differ- 
ential zero-frequency shot noise, as shown in Fig. 4. 
Differing from the differential conductance, the resulting 
QD spin-splitting causes an additional step structure in 
jpfc with equal (P) and unequal (AP) heights [Fig. 4(a), 
(c), and (d)]. Simple algebraic calculation gives: 



dS p _ dS AP 
v^o ~ dV 



dV 
dS p 



= 0, 



dV 
dS AP 



dV 



4 > 0, 
4(1- 



1 



> 0, 



(21) 
(22) 
(23) 



AP 



at A = 0. Generally, the zero- frequency shot noise will 
always increase with increasing bias-voltage. However, a 
contrary phenomenon is predicted in our calculation for 
the case of the AP alignment leads with large polarization 
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FIG. 5: Temperature dependencies of differential shot noise 
as functions of bias- voltage in the AP configuration without 
external magnetic field for p = 0.2 (a) and p — 0.8 (b). 
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FIG. 4: Differential zero-frequency shot noise as a function of 
bias- voltage with zero (a), (b) and nonzero (c), (d) external 
magnetic fields for various polarizations p. The other param- 
eters are the same as in Fig. 1. 



p: the differential shot noise changes sign near V = —A. 
This behavior can also be ascribed to the rather sharp 
peak in ^ r due to large p, i.e. to rapid QD spin-reversal. 
In Fig. 5, we also examine the temperature dependence 
of the differential shot noise at A = 0. It is seen that this 
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CONCLUSIONS 



In summary, we have applied a quantum Langevin 
equation approach to derive a Bloch-type equation micro- 
scopically for analytical studies of inelastic cotunneling, 
which we have carried out in detail, determining the lo- 
cal magnetization, current and its fluctuations, in a single 
QD attached to two ferromagnetic electrodes. Our stud- 
ies reveal a number of interesting novel characteristics 
intimately related to spin-flip processes in the AP config- 
uration: 1) a ZBP of the differential conductance at zero 
external magnetic field; 2) asymmetric peak-splitting in 
the presence of a nonzero magnetic field; and 3) decreas- 
ing shot noise near V = —A. In addition, we have found 
that these characteristics are robust against temperature. 
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